A quasi-local measure of inter-scale transfer: An approach to understanding 

turbulence 
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Many questions remain in turbulence research — and related fields — about the underlying physical 
processes that transfer scalar quantities, such as the kinetic energy, between different length scales. 
Measurement of an ensemble-averaged flux between scales has long been possible using a variety 
of techniques, but instantaneous, spatially-local realizations of the transfer have not. The ability 
to visualize scale-to-scale transfer as a field quantity is crucial for developing a clear picture of the 
physics underlying the transfer processes and the role played by flow structure. A general technique 
for obtaining these scale-to-scale transfer fields, called the filter approach, is described. The effects of 
different filters, finite system size, and limited resolution are explored for experimental and numerical 
data of two-dimensional turbulence. 

PACS numbers: 47.27-i 



I. INTRODUCTION 

Turbulence is characterized by the transfer of an invis- 
cid constant, such as kinetic energy or a passive scalar, 
between different length scales. One of the primary goals 
of turbulence research is to understand the mechanisms 
that drive this transfer process. Numerical simulations 
performed by Farge et al. suggest that in both two- 
dimensional (2D) and three-dimensional (3D) turbulence 
it is a relatively small number of "coherent structures" 
that dominate the turbulent dynamicsi^ As a result, 
it is important to understand how the existence and in- 
teraction of these coherent structures affect the transfer 
processes. Thus, the spatially-local transfer properties of 
the flow need to be measured and correlated with these 
structures. 

In this manuscript a tool for obtaining local infor- 
mation about the scale-to-scale transfer of inviscid con- 
stants from experimental/numerical data is examined. 
The method, called the "filter approach" (FA), has tra- 
ditionally been applied to large-eddy numerical simula- 
tions (LES) ^ but is developed here in the context of ex- 
perimental data analysis. By applying a low-pass spa- 
tial filter to the equations of motion (the incompressible 
Navier-Stokes equation), separate equations for the fil- 
tered, or large-scale, fields and the remaining small-scale 
fields can be derived^ Within the resulting equations are 
coupling terms that represent the interaction of the large- 
and small-scale fields with each other. In LES schemes 
the coupling terms in the large-scale equations are mod- 
eled thereby eliminating the necessity of directly com- 
puting the small scales. 3 In our application of the filter 
approach, however, data from direct numerical simula- 
tions (DNS) or high resolution experiments are used to 
directly evaluate the coupling terms and obtain a quasi- 
local measure of the inter-scale transfer 

For readers familiar with LES and turbulence model- 
ing it should be stressed that the philosophy driving our 



use of the filter approach is different from the a priori 
development of LES models discussed mi There the ob- 
jective was to determine empirically which of several LES 
modeling schemes for the large- to small-scale coupling 
terms most successfully emulates physical data. In that 
case, FA was used primarily as a benchmark, and only 
went as far as measuring the inter-scale coupling term. 
Rather than producing results for numerical benchmark- 
ing, FA can also be used as an analysis probe to deter- 
mine where and when in a flow scale-to-scale transfer of 
inviscid constants takes place. In this way one can isolate 
important interaction events and form an understanding 
of turbulence inter-scale transfer mechanisms. Of course, 
this understanding could eventually be incorporated into 
LES models. 

The filter approach is applied in this paper to both ex- 
perimental and numerical data. The intent is not to in- 
vestigate the underlying physics of the turbulence, which 
will be presented in later papers, but rather to determine 
the appropriate interpretation of the results and the lim- 
itations imposed by different filters, spatial boundaries, 
and finite measurement resolution. 



II. EXPERIMENTAL AND NUMERICAL 
SYSTEMS 

Experimental measurements were carried out in a flow- 
ing soap-film channel, a quasi-2D system in which decay- 
ing turbulence of low to moderate Reynolds number can 
be generated (10 2 < Re < 10 4 ). The channel was 5 cm 
wide and was inclined at an angle of 75° with respect 
to vertical. The mean flow was 120 cm/s and the film 
thickness was about 10 /im. A more detailed descrip- 
tion of the channel can be found in&S Using the em- 
pirical relationships measured in^S the films' kinematic 
viscosity was v w 0.03 cm 2 /s. The turbulence generat- 
ing grid consisted of rods of 0.12 cm diameter with 0.22 
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cm spacing between the rods. Thus, the blocking frac- 
tion is around 0.3, which is typical for turbulence in 2D 
soap film flows The resulting Reynolds number, 
Re = UL/v, was 880 based on the mean- flow velocity and 
an injection scale of Li n j = 0.22 cm. The turbulent veloc- 
ity, u(x), and vorticity, w(x), fields generated by the grid 
were obtained by tracking 3 — 5 (im polystyrene spheres 
(density approximately 1.05 g/cc) within a 1.8 x 1.8 cm 2 
region located 6 cm downstream from the grid (20-30 
eddy rotation times) The particles were illuminated 
with a double pulsed Nd:Yag laser and their images cap- 
tured by a 12-bit, 2048 x 2048 pixel camera. Around 
3 x 10 4 particles were individually tracked for each image 
pair and their velocities and local shears were interpo- 
lated to a discrete 135 x 135 grid. One-thousand ve- 
locity and vorticity fields were obtained in this way and 
were used to compute ensemble averages of the statistical 
measures described below. Typical velocity and vorticity 
fields are shown in Fig. ^ 

To supplement the experimental data, a direct numeri- 
cal simulation of the 2D Navier Stokes equation was per- 
formed. Computational details are presented in>i& The 
equation 
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was simulated in a square domain with side L = 2ir and 
periodic boundary conditions. Here, is the i-th. com- 
ponent of the velocity; uj — CijdjUi is the vorticity; and F 
is a stirring force applied to wave numbers |k| = 4 — » 7. 
The Einstein summation convention is used throughout. 
Two values of m are considered: m = 1 corresponding to 
Laplacian viscosity (v\ — 0.01), and m = 8 correspond- 
ing to hyper-viscosity {y% = 1.4 x 10~ 7 ), which has the 
effect of extending the inertial range. 

The equation was solved using a fully de-aliased, 
parallel pseudo-spectral code with second-order Adam- 
Bashforth time-stepping. The resolution was 2048 2 . A 
statistically-stationary state was achieved after about 
200 large-eddy turn-over times. Representative examples 
of vorticity fields generated using Laplacian and hyper- 
viscosity are shown in Fig. [21 



III. THE FILTER APPROACH 
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FIG. 1: Typical a) velocity and corresponding b) vorticity 
field obtained from the flowing soap film channel. The hash 
marks represent 1 mm increments, the mean flow (subtracted 
out) is in the — y direction. The top of the image is ~ 3 cm 
downstream from the energy injection grid. 



A. General Application 



The filter approach allows the direct measurement of 
the coupling between scales in non-linear systems. It was 
originally developed as a tool to truncate numerical sim- 
ulations of turbulence by modeling small-scale behavior 
with the knowledge of large-scale behavior, i.e. LES. 3 
In this paper, the general application of FA to nonlinear 
systems will be considered first, followed by a specific ex- 
ample: the case of energy and enstrophy transfer between 
length scales in 2D turbulence. 



General features of the transfer process are examined 
by considering a field, <?(x), that is evolved by a non- 
linear evolution operator, F^ nl \ as F^'q — 0. Given a 
scale, I, one can separate the field, q, into large-scale, qi, 
and small-scale, q s , components. This is done by convolv- 
ing a filter function, G/, with q and defining qi = G\ * q 
and q s = q — qi = (S — G{) * q, where S is the Dirac 
delta function. Applying Gi to the evolution equation 
for q yields Gi * (F^q) — 0. Adding zero, written as 
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F^qi - F^qi, results in: 

F^qi = -(Gi*(F (al) Q)-F (ia) Qi)- (2) 

Eq. J5J indicates that the large-scale field, qi, evolves in 
exactly the same manner as the full field, q, up to a cou- 
pling term C{q h q) = -{Gi * (F^q) - F^qi). This 
coupling term, which can be considered as an external 
forcing or damping, arises from the small-scale field, q s , 
interacting with q^ Thus, C provides a measure of the 
interaction between large and small scales. The mea- 
surement of C requires no assumptions about the field, 
q, such as homogeneity or isotropy. Moreover, C is a field 




FIG. 2: DNS vorticity fields with a) Laplacian viscosity (m = 
1) and b) hyper-viscosity (m = 8) 



quantity, not simply an average, and can reveal not only 
information about the coupling between scales, but also 
when and where such interactions are taking place and 
with what strength. 

The field nature of the coupling term derived above, 
G, makes FA very valuable in the context of studying 
turbulence. As mentioned in the introduction, turbu- 
lence is characterized by an average scale-to-scale flux 
of quantities such as the kinetic energy. Many possible 
mechanisms underlying the transfer processes have been 
suggested, e.g., the stretching of vortex tubes into thin 
filaments may account for some fraction of the down- 
scale transfer of energy in 3D, 1 ? In 2D, vortex merger has 
been postulated as a way of transferring energy to larger 
scales^ The difficulty is that neither of these pictures has 
been conclusively correlated with topological flow struc- 
tures or the underlying transfer dynamics, though some 
attempts have been made^ By using FA to measure the 
inter-scale transfer and correlating this field with the po- 
sition of flow structures (identified by other means), one 
can determine the veracity of these hypothesized transfer 
mechanisms. 

There are a number of subtleties to consider when 
applying FA to experimental data. First, there is the 
choice of the filter function, G;. The interpretation of 
the scale-to-scale transfer depends on the selection of the 
filter function. For example, if the convolution function 
is Gaussian defined by length scale I, the interpretation 
of the filtered functions is as given above. On the other 
hand, if the convolution is with the kernel Hi = (S — G/), 
where again G/ is Gaussian and S is a Dirac delta func- 
tion, the resulting interpretation of the convolved fields, 
Hi * /, as "large-scale" is incorrect. Indeed, the function 
Hi produces "small-scale" fields. Other ramifications of 
changing the filter function will be explored shortly but 
the convention is adopted that the filter function is al- 
ways low-pass. 

Two additional considerations when applying FA to 
real data are finite measurement resolution and the ex- 
istence of spatial boundaries. The extent to which the 
measurements are sensitive to either of these factors de- 
pends on the quality of the data and on the form of the 
filter used. Significant research in the context of LES 
has addressed similar concerns, but always for the pur- 
pose of approximating the physical system in a numerical 
simulation^ These issues are explored in some depth in 
a later section. 



B. Application to two-dimensional turbulence: 
energy and enstrophy transport 

In two-dimensional turbulence kinetic energy, E = 
u 2 /2, and enstrophy, £1 = uj 2 /2, are conserved in the 
inviscid limit. Theory, numerics, and experiments all in- 
dicate that energy is transferred on average from small to 
large length scales (up-scale) and that enstrophy is trans- 
ferred in the opposite direction (down-scale) FA, how- 
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ever, can be used to obtain much more detailed informa- 
tion about the energy and enstrophy transfer within the 
flow. We begin by examining the scale-to-scale coupling 
of energy in the 2D Euler equation (viscous terms are 
linear and, hence, cause no direct scale-to-scale transfer, 
eliminating the need to examine the full Navier-Stokes 
equation) : 



duj 
~dt 



duj 
'dxj 



dp 



(3) 



where Ui is the i-th component of the velocity field, p 
is the density normalized pressure field, and summation 
over repeated indices is assumed. Contracting this evo- 
lution equation with a filter function, Gi , and extracting 
the coupling term as described above yields: 



d(uj)i 
dt 



Mi 



d{uj 
Ox* 



dpi 



dx 



■((uiUj)i- (ui)i(uj)i). (4) 



The notation (/)/, or simply /), will be used to denote the 
large-scale field Gi * /. Eq. J3J is almost the equivalent 
of Eq. The one delicacy is that the term pi is not the 
large-scale pressure field, i.e., it is not the field obtained 
using the gradients of (ui)i. This is not an important 
issue here, however, because, as will be demonstrated 
later, pi does not contribute to the inter-scale transfer. 

From Eq. Q the scale-to-scale coupling term for the 
velocity is given by 



C((Ui) h Ui 



Or, 



dx; 



(5) 



where rjj — (uiUj)i — (ui)i(uj)i is the subgrid-scale stress 
tensor. To obtain the equation for the evolution of large- 
scale energy, — (ui)f/2 one multiplies Eq. Q by 
(ui)i. (Note that the notation is used rather than 
Ei since Ei = 
equation, 



)l/2 ^ (Ui)f/2 = £?(0.) The resulting 
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^ +^ 7 (K-)isW + (« J ), Pl ) = -(«i)j '" 
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dt 



dxj 



(6) 



for the energy contained at scales larger than I, is iden- 
tical to the full energy evolution equation up to the cou- 
pling term on the right hand side (again ignoring the 
pressure term). 

This coupling term is not yet in the form to directly 
yield scale-to-scale energy transfer information. There 
are two ways in which the small-scale velocities can 
change the large-scale energy: by physically transport- 
ing it from point-to-point or by transferring it between 
scales. To separate the two, the Leibnitz rule is used to 
rewrite the right hand side of Eq. J§} as 



_d_ 

dxj 



(') 



(7) 



Notice that the latter of these two terms is Galilean in- 
variant, whereas the former is not. Boosts to the refer- 
ence frame should not change the scale-to-scale transfer 
of energy but will change the point-to-point transport. 
Therefore, the former term is attributed to the point-to- 
point transport and the latter to the scale-to-scalc trans- 
fer. 

Another way to contrast the point-to-point coupling 
with scale-to-scale coupling is to consider the limit of 
a homogenous system. In this case taking an ensem- 
ble average, (...), should eliminate all point-to-point 
transport terms, leaving only inter-scale transfer contri- 
butions. Since the action of ensemble averaging com- 
mutes with the derivative operation, and since the spa- 
tial derivative of an ensemble average is zero in the limit 
of homogeneity, the ensemble average of Eq. © is sim- 
ply dt(E^) — (rjj' ] dj(ui)i) . Not only does this demon- 
strate that dj[r^ (ui)i] is a point-to-point term, but it 
also demonstrates that all of the terms on the left, other 
than the time derivative, are point-to-point as well. This 
fact allows us to ignore the delicacy with respect to the 
pressure term: it does not affect scale-to-scale transfer. 

For simplicity Eq. EJis rewritten as 



dxj 



where J^p and n^) are defined as 



n (0 



(i)d(uj)i 



v dx, 
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(8) 

(9) 
(10) 



where (Sij)i is the large-scale strain tensor. The negative 
sign in the definition of is added so that down-scale 
transfer has a positive value whereas up-scale transfer is 
negative. 

An almost identical method can be used to determine 
the scale-to-scale enstrophy transfer of the flow. The 
starting point, however, is the 2D Euler equation for vor- 
ticity, 



du 
~dt 



duo 
dxi 



0. 



(11) 



As above, the equation is contracted with a filter func- 
tion, G;, and the coupling term is extracted, 



duj, dun 

-dT + {u > )l dx- 



da 



(0 



dxj 



(12) 



where a^p — (ujco)i — (uj)iu>i is the subgrid scale vortic- 

ity transport vector. Notice that the term <j*p defined 

in Eq. ^|has an almost identical form to in the en- 
ergy equations. This general form is typical of the filter 
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FIG. 3: Obtaining the scale-to-scale enstrophy flux for a velocity field obtained from the soap film. The filter function, Gi, 
used was Gaussian with I = 0.2 cm. Row 1: Unfiltered velocity, u, vorticity, u), and vorticity transport ucu. Row 2: Filtered 
velocity, uj, vorticity, u>i, and vorticity transport, (uu>)i. Row 3: The subgrid vorticity transport vector a"', large scale vorticity 
gradient Vu>i and the scale-to-scale enstrophy transfer . 



approach. The coupling terms between large- and small- 
scale fields often take the form (ab)i — aibi for quadratic 
nonlinearities. 



To change the large-scale vorticity equation to an equa- 
tion for large-scale enstrophy evolution, one must multi- 



ply by uji. This yields 



dt 



+ 



dx, 



do 



(0 dui 
i dxj 



(13) 



where again the Leibniz rule was used to separate point- 
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to-point transport from inter-scale transfer. Grouping 
the appropriate terms, as was done for the energy equa- 
tion, leads to the final form, 



0fiW dK 



at 



(0 



dxj 



where 



Z« = 



_ (i) duii 
1 dx,j 



(0 



(14) 

(15) 
(16) 



where the negative sign in front of the scale-to-scale cou- 
pling term, Z"\ again makes down-scale transfer posi- 
tive. 

To make the above derivations more concrete, some 
examples are provided of filtered fields. For purposes of 
illustration, vorticity and enstrophy fields are presented 
and an analysis of velocity and energy fields is left to a 
later publication. Figure [3] displays steps in the calcula- 
tion of the scale-to-scale enstrophy transfer, Z^ l \ for a 
typical vorticity field extracted from experimental data. 
For this calculation the filter function was a Gaussian 
with Fourier-space definition 



G,(k) 



(17) 



where fc; = 2-k/I. These figures illustrate that the ap- 
plication of FA is straightforward: (1) compute the sec- 
ondary field, UiUi, from the measured velocity and vor- 
ticity fields; (2) perform a convolution of these fields to 



obtain erf 1 ; (3) take the scalar product of /j\ l> with the ap- 
propriate gradient of the large-scale fields, namely diU>i- 
There is, however, a caveat for the general case. It 
may not always be possible to separate the point-to-point 
transport from the scale-to-scale transfer terms. For the 
energy equation, the terms were determined by using the 
Leibniz rule to separate out the Galilean invariant part 
of the energy flux. For the enstrophy, the separation 
was obtained by analogy with the energy equation rather 
than by a strict application of Galilean invariance. There 
is no a priori expectation that such a separation will 
be as simple, or even possible, for arbitrary nonlinear 
systems. In compressible flows, for example, it is possible 
to measure the coupling terms, but the point-to-point 
transport caused by small scales is tied to the scale-to- 
scale transfer in a non-trivial way. The interpretation of 
such results must, therefore, be done carefully. 
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IV. SYSTEMATIC EFFECTS 

This section will address the effect of varying the fil- 
ter, Gi, on the interpretation of the results, the conse- 
quences of the data being limited in spatial extent {i.e. 
by boundaries), and variations in the results caused by 
experimental or numerical limitations on the resolution 
of the data. 




FIG. 4: Illustration of the various filter functions G, that 
are used in the paper (see Eq.EUand Eq.EHfor details) For 
both plots the real space convolution function is presented 
with the root of the Fourier spectrum inset, (a) Fourier filter 
convolution kernels, for orders n=2 (solid), n=3 (dash), n=4 
(dotted), n=5 (dash-dot), (b) Real space filters for orders 
n=2 (solid), n=3 (dash), n=4 (dotted) ,n=5 (dash-dot). 



A. Different Filters 

The interpretation of the results of FA may depend 
on the form of the filter function, G;. A Gaussian ker- 
nel was used in Eq. I|17|) . but this particular choice was 
made only because it has a simple interpretation in both 
real- and Fourier-space, taking the same form in both. 
The FA technique imposes no such constraints in general. 
In some cases, one may want to preferentially constrain 
the filter in real-space to a well-defined length scale or, 
instead, may want to filter so as to select only a sharp 
band of wave numbers. Either possibility can be explored 
using this tool. It must be kept in mind, however, that 
sharpening the filter in real- or Fourier-space causes a 
corresponding broadening of the filter in the other. The 
impact and interpretation of varying the form of the filter 
on the resulting fields is considered here. 

The form of the Fourier filters used is 

G| n) (k) =e"(^r)", (18) 
where n is the filter order. The case n = 2 is the Gaus- 
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FIG. 5: (a) Subsection of a typical vorticity field obtained from the soap film, (b)-(f) Calculation of from vorticity field 
shown in (a) using I = 0.15 cm and various Fourier and real space filters, (b) Gaussian filter (grayscale range ±1.95 x 10 7 s -3 ) 
(c) Fourier filter of order n = 3 (grayscale range ±2.47 x 10 7 s -3 ) (d) Real filter of order n = 3 (grayscale range ±1.93 x 10 7 s -3 ) 
(c) Fourier filter of order n = 4 (grayscale range ±2.83 x 10 7 s -3 ) (f) Real filter of order n = 4 (grayscale range ±1.90 x 10 7 s -3 ). 
The hatch marks represent 1 mm increments. 



sian filter considered earlier. As n increases, the filter 
sharpens around the Fourier mode corresponding to fil- 
ter length I. Similarly, the real-space filters are defined 
by: 

G ; (n) (x) = Ae-(w)". (19) 

Note that for n = 2, the real filter is equivalent to the 
Fourier filter. As the real filter is sharpened it approaches 
an area average over a box of diameter 21 /V. Both of 
these sets of filters are shown in Fig. Other types of 
filter are possible, but are not considered here. 

From an experimental point of view, the sharper real- 
space filters are more attractive than the Fourier fil- 
ters because the real-space envelope of the Fourier filters 
grows as the order of the filter is increased. Since exper- 
imental data is invariably windowed to the cross-section 
of the measurement apparatus, this means that sharper 
Fourier filters quickly grow to interact with boundaries. 
For sharp real filters this is not a problem: as the filter 
increases in sharpness it becomes more spatially com- 
pact. For the purpose of comparing with theory, however, 



sharper Fourier filters approach the ideal; see, for exam- 
ple, the discussion given in Frischii in which the fields 
are filtered by an infinitely sharp cutoff in Fourier-space. 

In Fig. [5] the results of using a Gaussian, two low-order 
Fourier-space filters, and two low-order real-space filters 
in the calculation of the enstrophy transfer are shown 
(the associated vorticity field is also displayed). Super- 
ficially, the fields are fairly similar. The strength of the 
fluctuations, however, changes with filter (the grey scale 
limits of the fields are noted in the captions) . The change 
in magnitude is stronger for Fourier filters than real-space 
filters, with the n = 4 Fourier-space filter experiencing 
50% larger fluctuations than for n = 2. 

The similarity in the fields seems to indicate that the 
relative magnitudes of enstrophy flux remain more or less 
constant. In particular, the stronger values of flux as- 
sociated with powerful vortices (middle top and middle 
left) have almost identical forms, though there may be a 
slight sharpening of the lobes for both higher order real- 
space and Fourier-space filters (which, we note, have a 
quadrupolar form) . Higher-order Fourier-space filters in- 
crease the symmetry of many features, i.e. they appear 
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less ovular (see for example the lobes indicated by the 
white arrows in Fig. [SJ, whereas the contrast between 
features and the background is enhanced by higher-order 
real-space filters. On the whole, however, the qualitative 
features are fairly insensitive to the form of filter used, 
maybe surprisingly so. One might expect that the beat- 
ing of the sharper Fourier kernel in real space would show 
up more strongly in the fields. For the stronger enstrophy 
transfer signals this does not appear to be the case. 

Weaker signals are more sensitive to the particular 
choice of filter. Consider the weak lobe in the lower mid- 
dle of the field indicated by the black arrow in Fig. 
(or the one in the middle right also indicated by a black 
arrow). For the Gaussian filter, the lobe is barely visi- 
ble. For sharper filters in both real- and Fourier-space 
(n = 4), however, the prominence of this lobe with re- 
spect to the stronger signals in the system is enhanced. 

The most important observation is that the qualitative 
structure of the field is relatively unaffected by the choice 
of filter. Sharper filters increase the "contrast" , but nei- 
ther eliminate nor create structural features in the flow. 
For the purposes of correlating scale-to-scale transfer to 
topological features, this is an important feature of FA. 

The above comparison is only qualitative. A more ac- 
curate quantitative comparison is presented in Fig. [ft] 
where the probability distribution function of the enstro- 
phy transfer is presented for the same set of filters used 
in Fig. \5\ The agreement between the different filters is 
quite good, but a little deceptive. Note, first, that the 
magnitude of the RMS fluctuations has been normalized 
out. Second, there may be (though it is, perhaps, below 
the noise level) a slight increase in asymmetry, in partic- 
ular in the negative tail of the PDF for real-space filters 
(the open symbols). Given that in Fig. the large values 
of enstrophy transfer, corresponding to the tails of the 
PDF, are relatively insensitive to changes in the filter, 
this collapse of the PDFs is reasonable. 

On the other hand, in Fig. the weaker values of 
inter-scale transfer have somewhat increased contrast. 
This change is emphasized in the lowest-order moments 
of the distribution, rather than in the tails. The aver- 
age, (Z^'), and fractional sign-probability comparison, 
6P(Z®) = P(ZM > 0) - P(Z® < 0), are shown in 
Fig. [7] Here, there is a significant difference in the av- 
erage enstrophy flux as a function of filter order. In 
particular, the average rises and falls more sharply for 
the higher-order Fourier-space filters than it does for the 
real-space filters. Also, the area fraction saturates at a 
smaller value, then falls more quickly. The signs of the 
average flux and the area difference, however, are quite 
robust, although the magnitudes seem to vary (in some 
places by a factor of two over this range of filter orders). 
The rise in the peak of average scale-to-scale transfer is 
reminiscent of ringing such as takes place in the Gibbs 
phenomenon. It is not clear whether or not this is the 
source of the change. 




FIG. 6: The probability distribution function (normalized 
by RMS) for Z (l) obtained using a Gaussian (solid line) two 
Fourier-space filters of order n = 3 (solid squares) and n = 4 
(solid circles) and two real-space filters of order n = 3(open 
squares) and n — 4 (open circles). The filter length was 
I = 0.2 mm. 



1. Interpretation 

At this point, one might ask: for which filter is the 
result closest to the "real" enstrophy flux? This ques- 
tion depends entirely on what one means by "enstrophy 
flux" . The scale-to-scale transfer between wave numbers 
is most closely approximated by higher-order Fourier fil- 
ters. Because of the associated broadening of the filter in 
real-space, however, the resulting fields are not as good 
a measure of the spatially-local enstrophy transfer. On 
the other hand, for the flux to be localized in physical 
space for comparison with real-space structures, sharp 
real-space filters arc preferable. In this limit, the enstro- 
phy flux can no longer be defined as the movement of 
enstrophy from small Fourier modes to larger ones but 
is actually a measure of the flux out of some bands and 
into others (which do not necessarily have larger wave 
number). The strength of FA does not lie in its ability to 
measure the exact magnitude of the transfer, but rather 
that the sign of the transport and the qualitative form of 
the fields is robust to changes in the filter. 



B. Boundaries 

An important consideration when using FA, one that 
is also a major issue in LES, is the presence of bound- 
aries (either physical or as limits of the viewing window) . 
As a Fourier-space filter grows in order, and correspond- 
ingly in real-space extent, less and less of the data can be 
used as the boundaries begin to affect the computation 
of the convolution in regions farther and farther into the 
interior. 

To investigate the influence of finite system size a ran- 
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FIG. 7: (a) The average enstrophy flux {Z (l) } and (b) SP(Z (l) ) 
(see text) for a range of length scales calculated using data 
from the soap film and a Gaussian filter function (solid line), 
two Fourier space filters of order n = 3(solid squares) and 
n = 4 (solid circles) and two real-space filters of order n — 3 
(open squares) and n = 4 (open circles). Inset in (a) is the 
average fluxes normalized by the RMS fluctuation size for the 
respective filters. 



dom periodic streamfunction was generated on a 512 x 
512 grid. This streamfunction was then used to obtain 
velocity and vorticity fields from which the enstrophy flux 
field, Z^p , was computed. The left half of the fields was 
then set to zero and the enstrophy field, recalcu- 
lated. The normalized RMS difference, 



Oh 



(^o^)rms 



(20) 



is shown in Fig. [5] as a function of x/l, the distance from 
the introduced boundary normalized by the filter size. 
The various plots are given for Fourier-space filters of 
different orders (see Eq. H5)l . 

With the exception of the order 3 filter there is a 
continuous increase in the boundary affects. Taking a 
nominal error rate of 10 -4 as acceptable, a boundary of 
b Ri l(n — 1) is appropriately sized for the data. This 



x/l 



FIG. 8: Interaction of various order Fourier-space filters with 
a boundary at x/l = 0. The graphs denote the RMS difference 
between the enstrophy flux measured in the presence of a 
boundary and that without. The various curves are for a 
Gaussian filter (solid- line), and Fourier filters of order n — 3 
(dashed), n = 4 (dash-dot), n = 6 (dotted), n = 8 (dash-dot- 
dot). 



results in a loss of 21 (n — 1) in linear box size since 
one must apply the condition to left-right (top-bottom) 
boundaries. This effective boundary has been adhered to 
throughout this paper. For real-space filters the bound- 
aries become less of an issue. Indeed, for the sharpest 
real-space filter (a simple average over a circle), the 
boundary is 1/2. 



C. Effects of Finite Resolution 

In this section the effect of finite measurement resolu- 
tion on the ability of FA to resolve the behavior of the 
enstrophy flux for a given filter function is explored. In 
the preceding sections experimental data was used in the 
analysis, but to estimate the effects of finite resolution it 
is necessary to consider numerical data where the range of 
scales and the uncertainty in the measured values is fairly 
well known. The numerical data also eliminates any con- 
cerns about the effect of boundaries since the boundary 
conditions are periodic. 

Pre-filtering at the limit of the data resolution should 
result in no apparent difference in the transfer from that 
computed from the original field. As the pre-filter length, 
L, is increased above the resolution scale differences in 
the measured transfer will increase. The question this 
raises is: how far above L are these finite-resolution ef- 
fects felt? To answer this question a series of 20 th or- 
der real-space filters with varying L are applied to two 
data sets simulating the local averaging inherent in ex- 
perimental measurements (such as image acquisition and 
particle tracking). From these pre-filtered fields the en- 
strophy transfer, Zjp , is calculated and compared with 
the full flux. 

The results of these calculations are shown for both 
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FIG. 9: The average enstrophy flux calculated using numer- 
ical data for (a) Laplacian viscosity and (b) hyper-viscosity 
with a pre-filter of length L = 0.01 (dashed), 0.02 (dotted), 
0.03 (dash-dot), 0.05 (dash-dot-dot). Also displayed is the 
total enstrophy flux (solid line) in the simulations. 

data sets in Fig. For both, the viscous cut-off scale is 
between 0.01 and 0.02. For pre-filters with L at or be- 
low this length there is little change in the magnitude of 
the average flux (slightly more prominent in the case of 
hyper- viscosity). Once the filter size exceeds the viscous 
scale, i.e., penetrates at all into the inertial range, the 
effects are quickly felt up to almost the injection scale, 
0.4. This may be a general feature of FA, or it may be 
a peculiarity of the 2D enstrophy transfer process being 
non-local. In either case, it demonstrates the importance 
of resolving the entire inertial range, including the vis- 
cous scale, for an accurate representation of the enstro- 
phy transfer to be possible. Whether or not enstrophy 
transfer is local in Fourier-space is an important question 
in its own right, and one which is sufficiently complex to 
warrant exploration in a separate paper. 

V. CONCLUSIONS 

In the above discussion and use of FA, familiar turbu- 
lence assumptions, such as homogeneity, isotropy, and in- 
ertial range, were not requirements of the measurement. 



This is perhaps the greatest strength of FA: it can be 
applied to any type of turbulent flow, regardless of that 
particular flow's properties. All that is needed is the 
evolution equation for the system from which the inter- 
scale transfer can be obtained. The flow doesn't have to 
be turbulent, FA can be applied to laminar or periodic 
flow and reasonable results obtained. In other words, the 
interpretation of FA does not rely on a pre-existing the- 
ory. Rather, it is a tool that can directly test notions of 
how inter-scale transfer takes place in systems and can 
be used to build appropriate theories. 

This "bottom up" approach to physics (rather than a 
theoretical trickle down) comes at a price: one must have 
highly resolved fields of data with a significant range of 
spatial scales. The former is a requirement imposed by 
the resolution issues discussed above, which can only be 
relaxed if one assumes a wavenumber local inter-scale 
transfer process. The latter is a limitation imposed by 
the data windowing and interaction of the filters with the 
measurement boundary. Of course, the standard tech- 
niques of measuring velocity fields in fluids, namely par- 
ticle imaging velocimetry or particle tracking, more or 
less ensure that the measurements are not far from the 
viscous scale, and thus of high enough resolution to not 
assume local transfer. This is because these techniques 
rely on groups of particles moving coherently, which only 
holds when the scales being probed are small enough 
that local Taylor expansions describe the flow. The dif- 
ficulty is in simultaneously measuring a significant range 
of scales above the viscous scale for meaningful informa- 
tion to be obtained. With standard PIV practices this is 
possible in 2D. Only with holographic PIV is this possible 
in 3D. There, however, the limitation is in the amount 
of data that can be obtained (order 10 fields is possible; 
currently, 1000 fields are not). 

As discussed in the introduction, our ultimate inten- 
tion for the FA technique is to probe the mechanisms 
driving the inter-scale transfer in 2D turbulence, with a 
particular eye towards the role of coherent structures in 
the transfer process. This will be done both via statis- 
tical analysis of fields in the laboratory frame of refer- 
ence (Eularian frame), as well as following the motion 
of fluid parcels (Lagrangian frame). FA, however, is not 
limited to these measurements alone, but could find ap- 
plication to the inter-scale transfer and mixing of passive 
scalars, or other quantities not necessarily conserved in 
the inviscid limit [i.e. local topology). And as hinted 
at previously, FA can be applied (carefully!!!) to general 
non-linear systems where an equation of motion is known 
and inter-scale transfer is of interest. 



VI. ACKNOWLEDGMENTS 

We thank Greg Eyink, Phil Marcus, Misha Chertkov 
and Boris Shraiman for interesting discussions and sug- 
gestions. 



11 



M. Farge, K. Schneider, and N. Kevlahan, Physics of Fluids 
11, 2187 (1999). 

M. Farge, G. Pellegrino, and K. Schneider, Physical Review 
Letters 87, 4501 (2001). 

M. Lesieur, Turbulence in Fluids (Kluwer Academic Pub- 
lishers, 1987). 

M. Germano, Journal of Fluid Mechanics 238, 325 (1992). 

B. Tao, J. Katz, and C. Meneveau, Journal of Fluid Me- 
chanics 457, 35 (2002). 

M. Rivera, W. Daniel, S. Chen, and R. Ecke, Physical 
Review Letters 90, 4502 (2003). 

C. Menevau and J. Katz, Annual Review of Fluid Mechan- 
ics pp. 1-32 (2000). 

M. Rivera, P. Vorobieff, and R. Ecke, Physical Review Let- 
ters 81, 1417 (1998). 

P. Vorobieff, M. Rivera, and R. Ecke, Physics of Fluids 11, 
2167 (1999). 

P. Vorobieff and R. Ecke, Physical Review E 60, 2953 



(1999). 

11 M. Gharib and P. Derango, Physica D 37, 406 (1989). 

12 H. Kellay, X. Wu, and W. Goldburg, Physical Review Let- 
ters 74, 3975 (1995). 

13 M. Rutgers, X. Wu, and W. Daniel, Review of Scientific 
Instruments 72, 3025 (2001). 

14 M. Ishikawa, Y. Murai, and F. Yamamoto, Measurement 
Science Technology 11, 677 (2000). 

15 K. Ohmi and H. Li, Experiments in Fluids 11, 603 (2000). 

16 S. CHcn, R. Ecke, G. Eyink, X. Wang, and Z. Xiao, Phys- 
ical Review Letters (2003). 

17 U. Frisch, Turbulence: The Legacy of A. N. Kolraogorov 
(Cambridge University Press, 1995). 

18 W. Daniel and M. Rutgers, Physical Review Letters 89, 
4502 (2002). 

19 R. Kraichnan and D. Montgomery, Reports on Progress in 
Physics 43, 547 (1980). 



